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ABSTRACT 

We present a large robust sample of 1503 reliable and unconfused 70 /im selected sources from the 
multiwavelength data set of the Cosmic Evolution Survey (COSMOS). Using the Spitzer IRAC and 
MIPS photometry, we estimate the total infrared luminosity, Lir (8-1000 /mi), by finding the best 
fit template from several different template libraries. The long wavelength 70 and 160 /im data allow 
us to obtain a reliable estimate of Lir, accurate to within 0.2 and 0.05 dex, respectively The 70 /im 
data point enables a significant improvement over the luminosity estimates possible with only a 24 /im 
detection. The full sample spans a wide range in Lir, Lir 10 s — 10 14 Lq, with a median luminosity 
of 10 11 Lq. We identify a total of 687 luminous, 303 ultraluminous, and 31 hyperluminous infrared 
galaxies (LIRGs, ULIRGs, and HyLIRGs) over the redshift range 0.01 < z < 3.5 with a median redshift 
of 0.5. Presented here are the full spectral energy distributions for each of the sources compiled from 
the extensive multiwavelength data set from the ultraviolet (UV) to the far-infrared (FIR) . A catalog 
of the general properties of the sample (including the photometry, redshifts, and Lir) are included 
with this paper. We find that the overall shape of the spectral energy distribution (SED) and trends 
with Lir (e.g., IR color temperatures and optical-IR ratios) are similar to what has been seen in 
studies of local objects, however, our large sample allows us to see the extreme spread in UV to near- 
infrared (NIR) colors spanning nearly three orders of magnitude. In addition, using SED fits we find 
possible evidence for a subset of cooler ultraluminous objects than observed locally. However, until 
direct observations at longer wavelengths are obtained, the peak of emission and the dust temperature 
cannot be well constrained. We use these SEDs, along with the deep radio and X-ray coverage of the 
field, to identify a large sample of candidate active galactic nuclei (AGN). We find that the fraction 
of AGN increases strongly with Lir, as it does in the local universe, and that nearly 70% of ULIRGs 
and all HyLIRGs likely host a powerful AGN. 

Subject headings: cosmology: observations — galaxies: active — galaxies: evolution — galaxies: high- 
redshift — infrared: galaxies — surveys 
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Luminous and ultraluminous infrared galaxies (LIRGs, 
L m = 10 u -10 12 Lq and ULIRGs, L m = 10 12 -10 13 L Q ) 
have played a major role in the study of galaxy forma- 
tion and evolution since their initial discovery. They 
were first identified in large numbers in the local uni- 
verse by the Infrared Astronomical Satellite (IRAS) and 
follow-up redshift surveys revealed that though they are 
rare in the local universe, their number density increases 
strongl y with redshift (e .g., the Bright Galaxy Survey 
(BGS): |Soifer et al.|1989| the Revis ed Bright Galaxy Sur- 
vey (RBGS): ISanders et al.||2003| Uy ULIRG Survey: 



Kim & Sanders 



rich of 



Since the launch of the Spitzer Space Telescope, many 
deep infrared surveys (particularly at 24 /*m where the 
MIPS detector is most sensitive) have been undertaken 
and have shaped our view of the significance of IR galax- 
ies out to z ~ 1. We now know that the cosmic star 
formation rate is dom inated by LIRGs beyond z > 0.7 
( Le Floc'h et al . 2005) and ULI RGs start to d ominate by 
z ~ 2 ( |Caputi et al. 12007] |Magnelli et al.|2009[ ). However, 
one limitation m all 24 /im selected studies is the diffi- 
culty in obtaining an accurate and reliable measurement 
of the total infrared luminosity of these sources since the 
24 /im selection wavelength represents rest-frame 12 /im 
by z — 1 and 8 /im by z — 2. Since the peak of the 
infrared emission for star-forming galaxies and galaxies 
with AGN is typically at ~ 50 — 200 /im, at higher red- 
shifts the selection wavelength moves further from the 
peak and becomes heavily influenced by the presence of 
polycyclic aromatic hydrocarbon (PAH) features. 

In order to study how various galaxy properties are re- 
lated to luminosity, an accurate estimate of th e total in- 



frared lum i nosity is needed. Previ ous studies (Bavouzet 
et al.|2008 Symeonidis et al.|2008[ ) have shown that rest- 



frame monochromatic luminosities in the mid-infrared 
(MIR) correlate very well with the total LrR, but they 
show a considerable amount of scatter. Of course, this 
problem becomes worse at higher redshifts as observed 
MIR probes shorter wavelengths. Ideally, one would use 
observations in the far-infrared (FIR) and submillimeter 
in order to sample the peak of the infrared emission in 
these sources, but current IR and submillimeter detec- 
tors lack the sensitivity to reach the necessary depths 
and only the most luminous sources are detected. Stud- 
ies at long wavelengths are also hampered by confusion 

due to the large beam si ze of most detectors . 

Previous studies (e.g., Huynh et al. 2007 Symeonidis 



et al.|2008||Magnelli et ai.|2009p based on 70 /7m selected 



sources have been hampered by this sensitivity problem 
and have lacked the depth or area necessary to detect a 
large sample of objects over a wide range of luminosities 
and redshifts. Deep surveys over small areas allow the 
detection of fainter, lower-luminosity sources, but they 
suffer from cosmic variance and poor statistics. Shal- 
low surveys over wide areas only sample the bright end 
of the luminosity function. In order to overcome these 
problems, we have obtained deep imaging at 24, 70, and 
160 fjan over the entire ~ 2 deg 2 of the COSMOS field 
(Le Floc'h et al. 2009, Frayer et al. 2009) . 
The Cosmic Evolution Survey (COSMOS: [Scoville 



et al.|2007b l was designed to study galaxy evolution over 



a large range of redshifts and environments, particularly 
with respect to large scale structure. The initial HST 



imaging (Koekemoer et al. 20071 has been complemented 
with wide multiwavelength coverage from space and the 
ground. Along with many spectroscopic redshifts, the 
deep multiwavelength photometry has allowed us to ob- 
tain photometric redshifts with unprecedente d accuracy, 
both for normal gal axies (Ilbert et al. 20091 and AGN 
( Salvato et al. 2009 1 . This unique data set provides the 



first opportunity to study the properties of a large sam- 
ple of 70 /im selected sources in detail over a wide range 
in redshift, luminosity, and environment. 

This paper is the first of a two part series on the prop- 
erties of a large sample of ^1500 70 /im selected galaxies 
in the COSMOS field. §2 summarizes the COSMOS data 
sets used for this study and our sample selection is de- 
scribed in §3 along with a summary of general properties 
of the sources. We present the full UV-FIR SEDs and 
describe the estimate of their total infrared luminosity 
in §4. We discuss our results and their implications in 
§5 and present our conclusions in §6. Paper II will dis- 
cuss the morphological and optical color properties of 
these sources as a function of redshift and luminosity. 
Throughout this paper we assume a ACDM cosmology 
with H = 70kms- 1 Mpc~ 1 , Q A = 0.7, and ft m = 0.3. 
All magnitudes are in the AB system unless otherwise 
stated. 

2. THE DATA SETS 

The imaging data used for this st udy were obtained a s 
part of the HST-COSMOS project rtScoville et al.|2007al 



COSMOS originated as an HST Treasury program imag- 
ing an ~ 2 deg 2 equatorial field with the Advanced Cam- 
era for Surveys (ACS), using the F814W filter (I-band). 
This is the largest contiguous field ever observed by HST. 
To supplement the ACS coverage, follow-up observations 
have been obtained across the entire spectrum using both 
space and ground based facilities. A complete descrip- 
tion of the C OSMOS data sets can be found in IScovillel 
et al. ( 2007b I . In this section we provide a brief descrip- 
tion of each data set and refer the reader to the relevant 
papers. These data sets are summarized in Table [I] 

2.1. S-COSMOS 



Spitzer-COSMOS (S-COSMOS: |Sanders et al.||2007 l 
is a Spitzer legacy survey designed to cover the entire 
2 deg 2 COSMOS field with both the MIPS and IRAC 
instruments. The deep IRAC data for the full COSMOS 
field was taken during cycle 2 (Jan 2006) with a total of 
166 hrs to cover the field in all 4 bands (3.6, 4.5, 5.6, and 
8 /im) down to a total depth of 1200 s. The 5 a sensitivity 
at 3.6 /im is 0.9 /iJy. For a detailed desc ription of the 
data r eduction of all S-COSMOS data, see |Sanders et al.| 
(20071. 

The MIPS coverage of the COSMOS field was taken 
during cycle 2 and cycle 3. The aim of the cycle 2 obser- 
vations was to obtain shallow (80, 40, and 8 s per pixel 
at 24, 70, and 160 /im, respectively) imaging over a large 
~ 4 deg 2 area centered on the COSMOS field. In addi- 
tion to the shallow imaging, deeper imaging was obtained 
for a small test region covering ~ 8% of the field. The to- 
tal depths obtained over this test region were 3200, 1560, 
and 320 s at 24, 70, and 160 /im, respectively, with a 5er 
sensitivity of 0.071, 7.5, and 70 mJy. Deep MIPS imaging 
was obtained for the entire field during cycle 3 (Jan-May 
2007 and Jan 2008) with a depth of 3400, 1350, and 273 s 
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at 24, 70, and 160 /im with a 5<r sensitivity of 0.08, 8.5. 
and 65 mJy, respectively. 
The cycle 2 and cycle 3 MIPS 24 /im data were com- 



bined and red uced using the MOPEX package ( Makovoz 
fc Khan|2005 1 to produce a final mosaic. The source cat- 
alog was produced using SExtractor ( |Bertin fc Arnouts 
1996J) for source detection and the DAOPHOT pack- 
age (|Stetso n|198 7 ) for measuring flux densities using the 
ad~F 



Point Spread Function (PSF) fitting technique. This fi- 
nal catalog is > 90% complete above 80 /^Jy and ~ 80% 
complete above 60 /iJy. Details on the 24 /jm obser- 
vations, data reduction, and source properties are pre- 
sented in Le Floc'h etal. (2009, in press). 

The raw 70 and 160 /xm data were reduced using the 
Germanium Reprocessing Tools (GeRT) and then com- 
bined using the MOPEX package. The source detection 
was done using the Astronomical Point-Source Extrac- 
tion (APEX) tools within MOPEX. The systematic pho- 
tometric errors are 5% and 12% at 70 and 160 ^m, re- 
spectively. Details on the 70 and 160 /jm observations, 
d ata reduction, catalo gs, and source counts are presented 



Frayer et al. (2009). 



2.2. Hubble Space Telescope images 

The HST/ACS images were taken during HST cycles 
12 and 13 as part of the COSMOS treasury survey. A 
total of 583 orbits were used to obtain a coverage of 
1.64 deg 2 in the F814W filter. The ACS images were 
taken between Oct 2003 and June 2005 and reach a depth 
of 27.2 mag (5a). Details of the ACS images, includ- 
ing their calibrati on and reduction, are given in |Koeke-| 
|moer et al.j ( |2007[ ). In addition to the ACS images, NIC- 
MOS parallel images were taken over ~ 6% of the ACS 
area in the F160W (H-band) filter down to a depth of 
25.9 mag (5 a). The high resolution of the ACS images 
(FWHM = 0.09") is necessary for the detailed morpho- 
logical analysis presented in Paper II. The NICMOS par- 
allel images (FWHM = 0.16") are particularly useful for 
studying the effects of morphological /c-corrections. Due 
to the requirement of morphological information for this 
study, we limited our sample to sources within the ACS 
area. 

2.3. Ground based UV/ Optical/ NIR Data 

The groundbased UV, optical, and near-infrared 
(NIR) coverage consists of data taken with Subaru- 
SuprimeCam (Bj,Vj,g + ,r + ,i + ,z + ), CFHT-Megacam 
(u*,i*), UKIRT-WFCAM (J), and CFHT-WIRCAM 
(Ks). The complete description of the data reduction of 
the optical and near-infrared (NIR) ground based data 
can be found in Capak etal. (2009, in preparation). 

The broadband optical data (Bj,Vj,g + , 
were taken in 2004-2005 using Suprime-Cam on the 8 m 
Subaru Telescope on Mauna Kca in Hawaii. The seeing 
ranged from 0.4-2.0"and the depth at i + is 26.2 mag. 
Deep u* and i* images were obtained using the Mega- 
cam camera at the 3.6 m Canada France Hawaii Tele- 
scope (CFHT) with a depth of 26.5 mag at u* . 

The J-band images were taken in several observing 
runs between 2004-2007 with the Wide-Field Camera on 
the 3.8 m United Kingdom Infrared Telescpe (UKIRT) 
and the Kg-band images were taken in queue mode be- 
tween 2005-2007 with the Wide-field Infrared Camera 



, r+, i 1 , z' 



4+ -?+\ 



(WIRCAM) at CFHT. The data were reduced using cus- 
tom scripts in . IRAF 21 . The TERAPIX software Scamp 
(Bcrtin 2006) was used to calculate the astrometric so- 
lution and the fina l images were combined using Swarp 
(iBertin et al. 2002). The final catalogs are complete to 
23.7 and K$ = 23.5 mag. For details on the re- 
duction of the K§ images, see McCracken et al. (2009, 
submitted to ApJ). 

Photometry for all of these data sets was produced 
using SExtractor in dual mode (IBertin & Arnouts 1996 1 



J 



using the I-band image as the detection image and have 
been combined into a master catalog. 

2.4. GALEX UV Data 

Ultraviolet data were taken using the Galaxy Evolution 
Explorer (GALEX) in 2004. Four pointings of 50 ks each 
were used t o observe the entire COSMOS field (Zamojski 
et al. 2007) in the near and far-ultraviolet bands (1NUV 
and FUV)! The u* images from CFHT were used as a 
prior for the flux measurements of the GALEX images. 
These data reach a limiting magnitude of 25.7 and 26.0 
in the FUV and NUV, respectively. 

2.5. X-ray Data 

The COSMOS field was observed with the XMM- 
Newton satellite in 2004-2005 for a total of - 1.5 Ms 



( |Hasinger et alTp007 |. A total of 1887 point sources were 
detected in at least one of the bands (soft: 0.5-2 keV, 
hard: 2-10 keV, and ultra-hard: 5-10 keV). The limiting 
fluxes in each of these bands is 5 x 10~ 16 , 2 x 10~ 15 , 
and 5 x 10~ 15 ergcm~ 2 s - 1 for the soft, hard, and ultra- 



hard bands, respectively (Cappelluti et al. 2009). The 
Chandra-COSMOS survey (C-COSMOS) observed the 
central 0.5 deg 2 of the COSMOS field to a depth of 160 ks 
and an outer 0.4 deg 2 t o a depth of 80 ks with the Chan- 
dra X-ray Observatory ( |Elvis et al. 2009). The flux lim- 
its obtained are 1.9 x 10 " and 7.3 x 10 -16 erg cm -2 s _1 
for the soft (0.5-2 keV) and hard (2-10 keV) bands, re- 
spectively. A total of 1761 point-sources h ave been de- 
tected and are presented in the catalog of |Elvis et al. 
(2009). In this paper, we use the optical identification 
catalogs (Brusa et al. 2009, Civano et al. 2009, in prepara- 
tion) from both of these data sets to identify AGN among 
our 70 /im selected sample of galaxies. 

2.6. Radio Data 

The COSMOS field was observed with the Very Large 
Array (VLA) at 1.4 GHz as part of a VLA Large Project 
in 2004-2005 for a total of 280 ho urs in A and C configu - 
ration to cover the entire 2 deg 2 ( |Schinnerer et al.|2007 |. 
The data has a mean rms noise ot 10.5 /iJy/ beam and a 
5 a sensitivity of 55 /iJy. The cat alog used here co nsists 
of - 2417 sources with S/N > 5 ( |Bondi et al.]|2008l ). 



2.7. Spectroscopy 

Spectra of the 70 /xm sources were obtained from sev- 
eral different sources. The first is a follow-up survey of 
24 /im sources in the deep MIPS test region of the field 
using DEIMOS on Keck II. The complete survey details 
will be described in Kartaltepe et al. (2009, in prepara- 
tion). The observations were taken on 5 nights between 

21 http://iraf.noao.edu/ 



4 



Kartaltepe et al. 



2007-2009 using the 600 1/mm grating with a 1" wide slit 
and typical slit lengths of ~ 6". The typical wavelength 
coverage with these settings is 4000-10000 A with a dis- 
persion of 0.65 A per pixel and a spectral resolution of 
4.6 A. The total integration times were between 1 and 3 
hours, depending on the magnitudes of the sources. The 
data were reduced using the DEEP2 pipeline 22 and red- 
shifts were measured individually for each source and a 
quality flag 1-4 was assigned, with 4 being a very secure 
redshift based on multiple lines, 3 being a likely redshift 
based on multiple lines, 2 being a less reliable redshift 
often based on only one line and a noisy spectrum, and 1 
being a best-guess redshift based on a noisy spectrum. A 
total of 62 70 /im detected sources have DEIMOS spec- 
troscopy. 

Spectra w ere also drawn fro m the zCOSMOS-bright 
10k sample (Lilly et al. 2007) observed with the Visi- 
ble M ulti-Object Spectrograph (VIMOS: |Le Fevre et al 
2003 I on the 8 m Very Large Telescope (VLT) in Paranal, 
Chile. This sample consists of the first 10,000 galaxies 
from a large flux limited (/ < 22.5) survey of galaxies 
over the entire COSMOS field. The zCOSMOS-bright 
survey used the R ~ 600 MR grism with a velocity reso- 
lution of < 100 km s -1 and a wavelength range of 5550- 
9650 A. 

Additional spectra have been obtained from a sur- 
vey of 677 XMM-Newton selected AGN candidates with 
the Inamori Magellan Areal Camera and Spectrograph 
(IMACS) on the 6.5 m Magellan (Baade) Telescope from 
2005-2008 on a total of 18 clear nights ( |Trump et ah 
2007||2009[ ). The 200 and 150 1/mm gratings were used. 
The wavelength range of these spectra is 5600 9200 A 
with a resolution of 10 A for the lower resolution 150 
1/mm grating. While AGN were the main targets of 
this survey, add itional sources were o bserved when pos- 
sible. Similarly, Prescott et al. (2006) conducted a spec- 
trocopic survey of quasar candidates using the Hectospec 
multiobject spectrograph on the 6.5 m Multi Mirror Tele- 
scope (MMT) on Mt. Hopkins using the 270 1/mm grating 
and covering a wavelength range of 3100-9000 A. Four 
sources were also observed using the Focal Reducer and 
low dispersion Spectrograph (FORS1) on the VLT (PI: 
C. Faure) covering a wavelength range of 3300-11000 A 
(Anguita et al. 2009, submitted). 

In addition to the surveys described above, spectra of 
the 70 /jm sources were found on the NASA Extragalactic 
Database (NED) p ublic archive from the Sloan Digital 
Sky Survey (SDSS: |Abazajian et al.||2009|) and the 2dF 
Galaxy Redshift Survey (2dFGRS: |Colless et al.|[200l l. 
See Table [2] for a summary of the spectroscopic surveys 
and a breakdown of the number of sources that were 
observed by each. 

3. THE 70 fj,M SAMPLE 

We restricted our study to the area of the field with 
HST/ACS coverage and selected sources at 70 /im with 
S/N > 3. After removing spurious sources detected 
around bright objects and removing sources in areas 
around bright stars masked out in the optical or IRAC 
photometry, the number of sources remaining is 1743. 
These 1743 sources represent the initial catalog of 70 /im 

22 http: //astro. berkeley.edu/~cooper/deep/spec2d/ 
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Fig. 1. — Histogram of 70 /im source fluxes in the deep MIPS test 
area (covering ~ 0.16 deg 2 : gray dashed) and the full COSMOS 
field (covering 1.64 deg 2 ) excluding the deep area (black). The test 
area represents ~ 9% of the COSMOS field and contains ~ 12% 
of the 70 fim sources (includes more sources since this region is 
deeper). The flux limit for the full field is ~ 6.5 mjy and the flux 
limit for the deep area is ~ 5 mjy for sources with S/N > 3. 

selected sources discussed below. This catalog includes 
sources within the deeper cycle 2 MIPS coverage area 
as well as the full field cycle 3 coverage, therefore, the 
deeper area has a slightly lower flux limit (~ 5mJy). 
The distribution of 70 /im fluxes for the deep area and 
the rest of the field is shown in Figure [l] 

3.1. Identification of Counterparts 

In order to ensure the most reliable identification of the 
optical counterpart for each source in our sample, we first 
matched the cleaned 70 /im catalog to the 24 /im catalog. 
Figure [2] shows the separation between the 70 /im source 
positions and the nearest (as well as the second and third 
nearest) 24 /im source. 79 of the 70 /im sources (4.5%) do 
not have a 24 /im counterpart within the 18" 70 /im beam 
or there is one on the very edge that does not appear to be 
a likely counterpart. The vast majority of these sources 
(92%) are at the low S/N end (3 < S/N < 4). In order 
to not be detected at 24 /im, these sources would have to 
have a 70 to 24 /im flux ratio of > 100. Sources with this 
extreme flux ratio are very rare in the rest of the sample 
with 24 /im counterparts (at the 1% level) . Physically, 
such a large flux ratio is possible at z > 2 where 24 /im 
could potentially pass through a silicate absorption fea- 
ture at rest frame 8/im, however, given the large number 
of sources not detected it is unlikely that all of them are 
such extreme high redshift sources. Furthermore, there 
are no sources in the local universe with such extreme 
ratios between rest-frame 8 and 24 /im. This, combined 
with their low S/N, indicates that they are most likely 
spurious detections and we therefore remove them from 
further analysis. 

Of the remaining 1664, 1528 sources were successfully 
matched to a 24 /im counterpart. Of the 1528 success- 
fully matched sources, 759 (~ 50%) have a single isolated 
24 /im source within the 70 /im beam making the identi- 
fication of the correct counterpart straightforward. For 
513 sources (~ 34%) there is more than one 24 /im source 
in the beam but the most likely counterpart is clear (i.e., 
the source closest to the center of the beam is also the 
brightest at 24 /im). 256 sources (~ 17%) have more than 
one 24 /im source in the beam but the likely counterpart 
is not as clear as in the previous case. For these sources 
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Fig. 2. — Top: Angular separation between each 70 /im source 
position and the nearest (as well as the second and third nearest) 
24 /im source. Bottom: Separation between matched source 70 
and 24 fim positions for the three different cases - those that are 
isolated (black), those that have a clear match (blue dashed), and 
those with a best-guess match (red dotted). 

we chose the source with the highest 24/iin flux as the 
best-guess match to the counterpart. Although it is pos- 
sible that some fraction of these "best-guess" matches 
are the wrong counterpart, we find that the results of 
our analysis do not change if they are left out. We esti- 
mate that our false match rate is on the order of 2.5%. 
For the 136 sources that cannot be matched there either 
appears to be more than one equally likely counterpart 
or the 70 /zm source appears to be a blend of multiple 
24 /im sources. It is highly probable that the multiple 
24 /im sources within the beam each contribute to the 
70 /im flux and therefore it is not possible to determine 
an accurate 70 /im flux for each source. Leaving these 
sources in our sample will only add an additional uncer- 
tainty, therefore we omit these ambiguous sources from 
further discussion. 

We then matched the remaining 1528 sources with be- 
lievable 24 /im counterparts to the IRAC source posi- 
tions. All but 9 of them have a clear match within 2" 
of the 24 /im source position. For these 9, the 24 /im 
source appears to be a blend of two IRAC sources so we 
remove them from the sample for the same reasons de- 
scribed above. For the remaining 1520 sources we then 
matched the IRAC counterparts to the optical I-band po- 
sitions and found clear counterparts for all but 11. Five 
of these appear to be two optical sources with blended 
IRAC photometry so we also exclude these 5 sources from 
our sample. The remaining six sources (0.3%) do not 
have an optical counterpart above the flux limit of the 
catalog. Without an optical counterpart we have no red- 
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Fig. 3. — Differential source counts (dN/dS X 5 2 5 ) for the fi- 
nal sample of 1503 sources compared to the published S-COSMOS 
(Frayer et al. 2009) and GOODS-N (Frayer et al. 2006) source 
counts. Note that the counts agree except at the lowest flux bin 
(570 < 6.5 /ijy). 

shift measurement for these six sources and thus cannot 
obtain an estimate of the total infrared luminosity so we 
exclude them from our final sample but note that their 
IRAC/MIPS properties do not appear to differ from the 
rest of the sample. Of the remaining 1508 sources, 5 are 
stars and so we remove them from the final sample as 
well. This results in a final sample size of 1503 70 /im 
selected sources with reliable counterparts at 24 /im, the 
IRAC bands, and the optical. We believe that the 9% 
of sources removed, mostly due to blending, will not af- 
fect our final result and removing them leaves us with a 
robust sample. 

We plot the differential source counts {dN/dS x S 2 ' 5 ) 
for the final sample of 1503 sources in Figure [3] along with 
the published source counts from S-COSMOs (Frayer et 
al. 2009) and GOODS-N (Frayer et al. 2006). The 
source counts for our sample match the published source 
counts very well except at the lowest flux bin (570 < 6.5 
jttJy). This indicates that spurious sources at low S/N do 
not have a significant contribution and that the sample 
is therefore reliable though incomplete at these low flux 
levels. 

The cross-matching at 160 /im was done by looking for 
160 /im sources near each of the 70 /im sources. For the 
vast majority of the cases (89%), only one 70 /im source 
was identified as a possible match to a 160 /im source. 
For most of the cases where more than one 70 /im source 
could have been identified as a possible match, one was 
much closer to the center of the beam and the other was 
far enough away that the choice was unambiguous (9%). 
For the handful of ambiguous cases (remaining 2%), we 
chose the brighter 70 micron source as the match. In 
total, we find that 463 of the 1503 70 /im selected sources 
(~ 31%) are detected at 160 /im. 

3.2. Photometric Redshifts 

All of the UV-NIR (including GALEX, IRAC, and 
all of the ground based optical and NIR) multiwave 
length photometry was used to pr oduce a photomet- 



ric redshift catalog for all galaxies dllbert et al.||2009 l 
and X-ray detected AGN ( jSalvato et al.||2009[ ) in the 
COSMOS field. Thanks to the depth and broad range 
of this multiwavelength coverage, the COSMOS pho- 
tometric redshifts have an unprecedented accuracy of 
dz/(l + z) = 0.007 at i+ < 22.5 and 0.012 at i+ < 24 
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Fig. 4. — Top: Redshift distribution of the entire 70 fim selected 
sample (white: photometic redshifts) and those with spectroscopy 
(black). Bottom: Fraction of sample with reliable spectroscopic 
redshifts as a function of redshift. Note that spectroscopic redshifts 
are available for ~ 40% of the sample and spread across the entire 
redshift range. 

out to z < 1.25. The AGN photometric redshifts have an 
accuracy of dz/(l + z) = 0.02 and 0.03 for i + < 22.5 and 
22.5 < i + < 25, respectively. The number of outliers for 
both galaxies and AGN are negligible. The accuracy and 
high precision of these photometric redshifts allows us to 
use our entire sample for the analysis of the properties 
of 70 fim sources instead of limiting our sample to those 
with spectroscopy as has generally been done in the past. 

The redshift distribution of the entire 70 fim sample is 
shown in Figure H] Photometric redshifts are available for 
all of the 1503 galaxies in the sample. 154 of these sources 
are detected in the X-ray by XMM-Newton and Chandra 
and so we use the photometric redshifts determined by 
Salvato et al.| ( 2009 ) . Our 70 /xm selected sample peaks at 
z ~ 0.35 and the number of sources drops beyond {z ~ 
1). However, given our large sample size, a significant 
number (273: 18%) are detected at 1 < z < 3 

3.3. Spectroscopic Redshifts 

Reliable spectroscopic redshifts are available for 602 
sources in the sample (~ 40%) and are represented by the 
filled histogram in Figure [4] Golumn 5 in Table [2]lists the 
number of sources in our sample that have spectra from 
each of the different redshift surveys. Several sources 
have redshifts available from more than one source. For 
these objects we compared each available redshift to 
check their reliability and chose the best one. A compari- 
son between the photometric and spectroscopic redshifts 
for those sources with very secure redshifts (quality flags 
3 and 4) is shown in Figure [5] The XMM sources are 
shown a s the red points and use the photometric red- 



shifts of Salvato et al. ( 2009[ ). The overall dispersion in 
this infrared selected sample is 0.02 x (1 + z). There are 
only 9 catastrophic failures (|z p — z s \/(l + z s ) > 0.15) in 
the sample (1.7%). For the rest of the paper, we adopt 
the spectroscopic redshift for those sources with reliable 
spectroscopy and the photometric redshift for the rest of 
the sample. 

3.4. UV-to-FIR Spectral Energy Distributions 

Figure [6] shows a compilation of all 1503 spectral en- 
ergy distributions ranging from the UV through the FIR 
binned by redshift and normalized at 1 /im. Over plot- 
ted in blue is the median SED for each luminosity bin. 
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Fig. 5. — Comparison of photometric and spectroscopic redshifts. 
Photometric redshifts are from llbert et al. 2009 (for normal galax- 
ies) and Salvato et al. 2009 (for X-ray detected AGN marked with 
red boxes). Note that the agreement is excellent for the AGN as 
well as the normal galaxies. The overall dispersion is 0.02 X (1 + z) 
(shown as dashed lines in the bottom panel) and the fraction of 
catastrophic failures is 1.7%. 

All of the photometry used to construct these SEDs, the 
redshifts, and the quantities derived in this paper (such 
as Lir) are listed in Table [3j 

3.5. X-ray Properties 
Using the XMM- Newton and Chandra X-ray catalogs 



( |Brusa et al.| |2009| Civano et al. 2009, in preparation) 
of the COSMOlTiield we identified which of the 70 /an 
selected sources are detected at X-ray energies. 119 of 
the sources were detected by XMM and 103 by Chan- 
dra with 68 detected by both, for a total of 154 X-ray 
detected sources. The spatial distribution of the sources 
detected by both XMM and Chandra is shown in Fig- 
ure [7] This plot illustrates the relative areas observed 
withboth telescopes within the ACS coverage area. The 
rest-frame X-ray luminosities for each of these sources 
was calculated from the hard band (0.5-10 keV) fluxes 
with an assumed power law photon index of T = 1.7 and 
is shown in the top panel of Figure [8] as a function of red- 
shift. 132 (85%) of these sources have Lx > 10 42 ergss _1 
(dashed line in Fig. |8| and so are likely to be AGN. 

Figure [9] shows the distribution of hardness ratios 
{{H - S)J[H + S), where H = 2 - lOkeV and S = 
0.5 — 2keV) for the 102 X-ray sources in the sample that 
are detected in both the soft and hard X-ray bands. For 
the remaining objects, 24 arc completely hard (i.e., no 
detection in the soft band), 29 are completely soft (no 
hard band detection) and only upper limits are available 
for 3. Over plotted on this histogram (dotted line) is 
the normalized distribution for the entire X-ray selected 
population that is detected in both bands. These two 
populations show striking differences. Whereas the full 
X-ray population has a hardness ratio that peaks around 
—0.5, the 70 /an selected population appears bimodal, 
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Fig. 6. — Spectral energy distributions of all 1503 70 fira sources binned by redshift and normalized at 1 fim. The SEDs are color coded 
such that the ones with the bluest slope in the UV-optical are purple/blue and the reddest are red and the spectrum in between shows the 
range of SED colors. The median SEE) for each bin is overplotted in blue. Note that at higher rcdshifts the SEDs extend further into the 
UV while the FIR points probe shorter wavelengths 
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Fig. 7. — Spatial distribution of all of the X-ray sources detected 
by XMM-Newton (dots) and Chandra (boxes). Note that a smaller 
area of the field was observed by Chandra. The two rectangles illus- 
trate the different depths of Chandra coverage. The outer rectangle 
(solid box) is ~ 0.9 deg 2 (covering 55% of the ACS area) repre- 
senting the entire area observed by Chandra. The inner rectangle 
(dotted box) is ~ 0.5 deg 2 (covering 30% of the ACS area) ob- 
served to an effective exposure time of ~ 160 ks. The area between 
the two rectangles was observed to an effective exposure time of 
~ 80 ks. 
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Fig. 8. — Top: Distribution of the rest-frame X-ray luminosi- 
ties for the 154 sources detected by either Chandra or XMM as 
a function of redshift. The dashed horizontal line represents the 
luminosity limit above which objects are considered to be AGN 
(Lx > 10 42 ergs s —1 ). Bottom: Distribution of rest-frame 1.4 GHz 
radio luminosities for the 562 sources detected by the VLA as a 
function of redshift. 
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Fig. 9. — Normalized distribution of hardness ratios for 70 /im 
selected sources (solid line) and the general X-ray population (dot- 
ted line) . Note that the 70 (im sample matches the general X-ray 
population very well at low hardness ratios but that there is an 
excess of 70 /im sources with larger hardness ratios. 

with peaks at ~ —0.5 and 0.35. The 70 fim. selected pop- 
ulation of X-ray detected AGN is on average harder than 
the general X-ray selected population. The large hard- 
ness ratios of some of the 70 sources indicates that 
there is significant obscuration present in these galaxies 
absorbing UV and soft X-ray emission. 

3.6. Radio Properties 

562 of the 70 /Ltm sources were matched to a radio 
source within 1 arcsecond. For each of these sources, we 
estimated the rest f rame 1.4 GH z luminosity (Li.4GHz) 



using equation 5 of Yang et al. (20071 and applying a 
fc-correction assuming a power law of the form F v oc u~ a 
where a, the ra dio spectral index, is assumed to be ~ 0.7 



( Condon|1992 1 . The radio luminosity is presented in the 
bottom panel of Figure [8] as a function of redshift for all 
562 radio detected sources in our sample. 

Using this subsample of objects detected in both the 
radio and infrared, we investigated the infrared-radio cor- 
relation. Even though this correlation has been shown to 
hold over a wide range in redshift and luminosity (e.g., 
Helou et al.||19851 ICondon et al.||1991a||Yun et ai.||2001| 
Appleton et al.||2004| |lbar et al.||2008p , there is some 
evidence that it might break down at extreme infrared 
luminosities or at hi gh redshift (e.g., i Kovacs et al.||2006 
Sajina et aL]|2008| (Younger et at] |2UU9| ) . A complete 
analysis of the evolution of the infrared-radio correlation 
with redshift is beyond the scope of this paper, how- 
ever, a full treatment of this issue will b e discussed in 
Sargent etal. (2009, submitted to ApJS). Condon et al. 
( 1991b I found the mean value of q for BGS starbursts to 
be < q >= 2.34 with a very small dispersion (a ~ .19) . 



This i s in ex cellent agreement with the results from 
et al.| ( |2001[ ) who found < q >= 2.34 for a sampli 
1809 g alaxies from the IRAS 2 Jy sample. The |Yun et al. 
(2001 1 study also found that there were slight variations 



in q with IR luminosity. For faint IR sources, they found 
values of q that were slightly higher than the mean and 
for bright IR sources they found a higher overall disper- 
sion (a ~ 0.33). 
We computed the value of q as defined originally by 



Helou et al.| ([1985J as the ratio between the FIR (a 
combination of infrared flux measurements from 60 and 
100 /im) and the flux at 1.4 GHz. We adopted this def- 
inition to allow for a direct comparison with previous 
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Fig. 10. — Radio/infrared flux ratio (g) as a function of total in- 
frared luminosity (top) and redshift (bottom). The horizontal line 
represents the mean value and the dotted lines are the values ±1 a 
from the mean. The scatter increases with LrR and the distribu- 
tion appears to bifurcate at ~ Ljr > 10 12 Lq . The large scatter 
at high luminosities is observed even when a more stringent (i.e., 
S/N > 5) cut is used. 



work. For our sample of 562 radio detected 70 /im se- 
lected sources we compute q using the rest frame 60 
and 100 /im flux densities obtained from the best-fit SED 
template and the fc-corrected 1.4 GHz flux density. The 
res ultin g q distribution for our sample is shown in Fig- 
ure 10 as a function of Lir (top) and redshift (bottom). 
We find a wide range in q (0.65-3.7) with an average 
q of 2.36 and a dispersion of 0.39 overall, in agreement 
with previous studies. However, we find that this disper- 
sion increases dramatically with Lir with an increase 
in both radio-excess (small values of q) and infrared- 
excess (large values of q) galaxies. For the LIRGs and 



ULIRGs in the sample, we find < q >= 2. 34 ±0.31 and 
< q >= 2.42 ± 0.50, respectively. It is thought that the 
slightly higher mean q at high luminositi es is the result 
of incr edibly high radio free-free opacity (Condon et al. 
1991b I. However, we note that this effect should decrease 
with -Lir in our sample as the rest frame radio frequency 
probed increases, and therefore, the infrared excess ob- 
served in these sources is not likely to be caused by this 
opacity effect. 

4. TOTAL INFRARED LUMINOSITIES 

4.1. SED Fits 

We obtained the total infrared luminosity, Ltr (8— 
1000 /im), for each source in our sample by fitting the 
8, 24, 70, and 160 /im data points to various SED tem- 



plates using the SED fitting code Le Phare 23 written by 
S. Arnouts and O. Ilbert. We weighted each of the four 
points equally to avoid over- weighting the points with the 
smallest error bars, which are typically those at shorter 
wavelengths (i.e., 8 and 24 /mi). The 8/im data point 
was only used at low redshifts. Once it shifted below 
7/im, it was no longer used due to confusion with PAH 
features and since the 24, 70, and 160 /im data points 
well represent the portion of the SED that contributes 
the bulk of the infrared luminosity. Beyond z ~ 2.4 the 
24 /im point is no longer used for the same reason. 

4.2. Template Libraries 
We us ed four different temp l ate libraries for the SEP 



fitting - Chary & Elbaz 


(2001), 


Dale & Hclou 2002, La- 


gache et al. 


(2003), and 


Siebenmorgen &; Kriigel| rt2007 ) 



resent a wide range of SED shapes and luminosities and 
they are widely used in the literature. Here, we briefly 
describe how each of these libraries was derived and dis- 
cuss the differences between them. For more information 
on these template librarie s, see the original paper s and 



t he detailed descr i ption in |Sym eonidis et al. (|2008 
: TEi 



is template library consists 



apt 

Chary & Elba^ \2001]) : 
of 105 templates based on the SEDs of four "prototypi- 
cal" galaxies (Arp220: ULIRG, NGC 6090: LIRG, M82: 
starburst, and M51: n ormal star form ing galaxy). They 
were derived using the |Silva et al.| ( |T998[ ) models with the 
MIR region replaced with ISOCAlvl observations. These 
templates were then divided into two portions (4-20 fjxa 
and 20-1000 /im) and interpolated between to arrive at 
a set of libraries o f varying shapes and luminosities. The 
Dale et al. ( 2001[ ) templates are also included in this set 
t o extencTthe range of s hapes . 

Dale & Helou\ ]2002\): These templates are updated 
versions of the |L)ale et ah J 2001 ) templates which used 
Desert et al.| ( |1990| ). This model in- 
large dust grains in ther- 
mal equilibrium, small grains semistochastically heated, 
and stochastically heated PAHs. They are based on 
IRAS /ISO observations of 69 normal star -forming galax- 
ies in the wavelength range 3-100 /im. |Dale fc Heloul 
(|2002|) improved upon these models at longer wave 



the model of Desert et a 
volves three components, 



length s usin g SCUBA observations of 114 BGS QSoifer 



et al. 19891 galaxies, 228 galaxies with ISOLWS (52- 



170 /mi; Brauher 2002) and 170 /im observations for 115 



galaxi es from the ISOPHOT Serendipity Survey ( |Stickel 



et al. 20001. All together, these 64 templates span the 
IR lumino sity range 10 s — \Q 12 L & 



Lagache et al 
built 



plates was 

uses the dust emission model of 
with a few modifications applied 



(2003): This library of starburst tern 
on the model of Maffei (1994) , which 



Desert et 
First 



al.| p90"l l, 

they replace 



the PAH portion of the templates with the one from 
Dale et al. (20011, then they increase the proportions of 
PAHs and very small grains while adding in extinction, 
and finally, they broaden the FIR peak and flatten the 
SED at long wavelengths. This library also includes tem- 
plates of normal star-forming galaxies derived from the 

the 



2000), 



ISOPHOT serendipity survey (Stickel et al. 
FIRBAC K galaxy SEDs, and long wavelength data"taken 
fromlDunne et al.l (120001) andlDunne & Ealesl (120011). 



http: / /www.cfht.hawaii.edu/~arnouts/LEPHARE/cfht_lcphare/lephare.h 
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Fig. 11. — Spectral energy distributions for a sample of 14 sources at a variety of luminosities and redshifts (both photometric and 
spectroscopic). Over plotted are the best fit templates from each of the four infrared template libraries as well as the best fit to the stellar 
component using empirical SEDs from Coleman et al. (1980). In order to illustrate the range of fits, half of the sources here have 160 /im 
detections and the other half do not. Note that at high redshift both the 70 and 160 fim points are on the low wavelength side of the peak 
and that the fits get progressively worse. Also, the fit to the stellar component becomes worse for sources with a significant contribution 
from AGN activity in the NIR (e.g., sources 1237 and 1539). 



: This template library consists of over 7000 model 
SEDs derive d using a radiat ive transfer and dust model 
described in Kruegel (2003). This model uses the as- 
sumption of spherical symmetry and divides the star- 
burst luminosity into two classes, that from OB stars in 
dense clouds (Lob) and the total from all other stars 
(Ltot — ^ob)- Five parameters are varied in the calcula- 
tions to create the grid of template SEDs: 1) The total 
IR luminosity ranges from 10 10 — 10 14 L Q , 2) The radius 
of the starburst nucleus: 0.35, 1, 3, 9, and 15 kpc, 3) 
The visual extinction of the nucleus: Ay =2.2, 4.5, 7, 



9, 18, 35, 70, and 120 mag, 4) The ratio of the luminos- 
ity contributed by the OB stars to the total luminosity: 
Lon/Ltot = 0.4, 0.6, and 0.9, and 5) The dust density in 
the hot spots: n hs = 10 2 , 10 3 , and 10 4 cm -3 . Combina- 
tions of parameters that are considered to be physically 
unlikely arc left out of the final grid. 



4.3. L 



in 



For each of the above template libraries, the best fit 
model was chosen by finding the one with the lowest x 2 
value and allowing for rescaling of the templates. The 
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Redshift (z) 

Fig. 12. — Total infrared luminosity (Ltr) as a function of red- 
shift for all 1503 70 /im sources in the COSMOS field determined 
from the best fit template using 8, 24, 70, and 160 fim where avail- 
able. The error bars represent the 1 a uncertainty for each source 
based on the fit, the redshift errors, and a 0.2 dex systematic un- 
certainty for the sources without a 160 (im detection. 

total infrared luminosity was then calculated from the 
best fit template by integrating from 8-1000 /im. A sam- 
ple of SEDs and the best fit from each template library 
is shown in Figure ITT] for a range of luminosities. We 
visually inspected the fits and found that they typically 
work just as well for AGN as for non-AGN (in the 8- 
1000 /im range) since our fit does not include all of the 
IRAC data points which would be the most affected by 
an AGN. For the final Ltr we adopt the value of the best 
fit of all t he template libraries. This final value is shown 
in Figure [12] as a function of redshift and is presented in 
Table [3] Since the template library is so large, a range of 
templates can be used to adequately fit the data points. 
We use this range, along with the 1 a uncertainties in 
the redshift measurement, to calculate the uncertainty 
on the final derived Ltr given in Table [3] In total we 
identify 687 LIRGs, 303 ULIRGs, and 31 HyLIRGs. 

4.4. Effect of 160 [im Detections on Ltr 

The depth and breadth of the 70 /zm imaging in the 
COSMOS field allow us to obtain a more accurate mea- 
surement of the total infrared luminosity than has pre- 
viously been possible for samples selected at 24 /im. The 
long wavelength selection means we can be confident in 
our estimate even at high redshifts since we are still prob- 
ing the infrared slope of the SED (35 /im at z = 1 and 
23 /im at z = 2). In addition, this large sample size of 
sources, including many that are also detected at 160 /im, 
enables us to statistically quantify how well we could 
have measured the luminosity if we only had detections 
in a subset of the infrared bands. This information is par- 
ticularly useful for comparing our results to other studies 
of sources selected at 24 /im. 

Figure [13] shows a comparison of our estimate of Ltr, for 
the subsample of sources detected at 160 /im (463 sources 
in total) with the estimate excluding the 160 /im data 
point (based only on 8, 24, and 70 /im). The solid line 
represents a 1-to-l relationship and the lower panel shows 
the difference between these two measurements. Since 
the majority of 70 /im selected sources are not detected at 
160 /im, it is important to understand and quantify this 
relationship. While the internal scatter of the points is 
quite small (a = 0.20 dex), there appears to be a signifi- 
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Fig. 13. — Comparison of the total infrared luminosity for the 
463 sources detected at 160 /im obtained with and without the 
160 /im data point. The error bars are the 1 a uncertainty based 
on the template fit. The bottom panel shows the difference in Lyr. 
between these two measurements. The mean offset is shown as the 
dotted line and the dispersion about this mean is given. 

cant offset (of 0.20 dex on average) between the estimates 
with and without the 160 /im data point, which increases 
with Ltr,. This offset indicates that the Lir determined 
without the 160 /im data point is underestimated by 0.20 
dex on average. 
In an attempt to isolate the possible cause of this off- 



set, Figures uA\ and [15 



contain the same plot divided 
into several flux and redshift bins, respectively. The off- 
set is largest for the sources at the faintest flux levels 
(0.36 dex at 4 < L(70) < 8mJy) and is negligible at the 
brightest fluxes (0.11 dex at F(70) > 20mJy). When 
divided into redshift bins, the largest offsets are seen at 
the high redshift end. The mean offset and internal dis- 
persion for each redshift and flux bin are summarized in 
Table [4] There are a few possible explanations for the 
observed offset. One is that the intrinsic noise in the 
160 /im measurement at these low flux levels biases the 
fit toward higher estimates than if the source was not de- 
tected at 160 /im at all. However, this explanation seems 
unlikely since the offset remains even when using higher 
S/N cuts at both 70 and 160 /im. Another possibility is 
that at higher redshifts the 70 /im detection samples too 
short a wavelength to put much of a constraint on the 
best fit template. This would result in an overall larger 
scatter but not an offset. And finally, it is possible that 
the detection at 160 /im preferentially picks out sources 
with a colder dust component than the rest of the 70 /im 
population. 

To test our fitting routine we fit the same data points 
(8-160 fim from IRAC and MIPS) for a set of 67 local 
LIRGs and ULIRGs fro m the Great Observa tory All-Sky 
LIRG Survey (GOALS: |Armus et alT1[2009| . As for the 
COSMOS sources, we obtained the best fit Lir with and 
without the 160 /im point and plot the comparison in Fig- 
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Fig. 14. — Same as Figure [T3] but divided into four 70 /im flux bins. Note that the mean offset between the two luminosity estimates 
decreases at higher 70 (im fluxes while the overall scatter stays about the same. 
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Fig. 15. — Same as Figure [13] but divided into seven redshift bins. Note that the mean offset between the two luminosity estimate 
increases with redshift while trie overall scatter stays about the same. 
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Fig. 16. — Comparison of the total infrared luminosity obtained 
with and without the 160 /im data point for a sample of 67 lo- 
cal (U)LIRGs from the GOALS survey. The low dispersion ob- 
tained between these two measures indicates that the template 
fitting method is reliable and that the rest-frame 70 fira data point 
is long enough to obtain a reliable measure of Ltr at low redshift 
(z < 0.1). 

ure [16] For these local sources we do not see the same 
offset that we see in the COSMOS sources and the intrin- 
sic scatter is much smaller (a — 0.08 dex). In addition, 
the Lir obtained from this fit agrees very well with Ltr 
in the literature from the IRAS data points. This test 
gives us confidence that the offset we observe is not due 
to a problem with the fitting routine and also indicates 
that the rest-frame 70 fim data point is at a long enough 
wavelength to obtain an accurate measurement of Ltr 
at low redshift. We also find that the template libraries 
based empirically on local sources tend to fit the data 
points better, as is expected. As an additional test we 
fit only the 8 and 24 /im data points for the local sources 
(roughly corresponding to observed 24 and 70 at z ~ 2) 
and find an overall scatter of ~ 0.17 dex. This result is 
consistent with the typical uncertainty we obtain for the 
COSMOS sources. 



pymeonidis et al. ( 2008 1 performed a similar compar- 
ison with 43 sources from the Extended Groth Strip. 
They derived an em pirical relation (Equation 1 in |Syme^| 
onidis et al. 2008) to obtain Ljr directly from the 24, 
70, and ISO /im tlux measurements without the need to 
fit templates or apply fc-corrections. This relation was 
found by summing the approximate contributions to the 
total infrared luminosity obtained using triangles and 
polygons of heights determined by the fluxes at the three 
MIPS wavelengths. They then add in a correction for the 
systematic offset found between this approximation and 
their total Ltr found through template fitting to arrive at 
the final equation. They find that this equation agrees 
well with their 43 EGS sources at all Ltr with a very 



Fig. 17. — Comparison of the total infrared luminosity obtained 
with the fits to the four infrared data points to that obtained with 
equation 1 of Symeonidis ct al. 2008. The agreement between our 
estimates and the equation is excellent although the dispersion 
increases at Lir > 10 12 Lq . At these higher luminosities, and 
consequently higher redshifts, the 160 fim data point is not long 
enough to precisely constrain the peak and thus there is a larger 
spread among the best fit templates. 



small scatter (0.06 dex). Figure 17 shows a comparison 
of the Ltr der ived here for the COSMO S sources and 
that using the Symeonidis et al. ( |2008[ ) equation. We 
also find excellent agreement and a small scatter (0.06 
dex overall), however, we find that the scatter increases 
beyond 10 12 L Q (~ 0.1 dex). A closer inspection of the 
sources with the largest deviation from the equation in- 
dicates that at these redshifts, the 160 /xm point is not 
long enough to constrain the peak of the emission and 
therefore, a wider range of templates fit the points than 
at lower redshifts. 
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Not detected at 160 fim _ 

Detected ot 160 ^tm . 













70 M m F" 1 ^* [mJy] 

Fig. 18. — 70 /im flux distribution of the 463 sources detected at 
160 fim (gray dashed) and the 1042 not detected at 160 fim. On 
average, those detected at 160 fim have higher 70 ^m fluxes. 
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Fig. 19. — Same as Figure [14] but with the luminosities obtained by fitting the mean SED for the 1042 sources not detected at 160 /im 
with stacked 160 (im flux measurements divided into several flux and redshift bins (shown as colored points). Note that the while the 
stacked points show the same amount of scatter overall between the two luminosity estimates, they do not have the same systematic offset 
as the 160 fim detected sources. This suggests that there is an intrinsic difference in the shape of the SEDs between sources with a 160 fim 
detection and those without. 



4.5. Sources not detected at 160 fim 

For the remaining 1042 sources without a detection at 
160 /im, the relationship shown in Figures |13f|15| indi- 
cates that an empirical correction may be needed. How- 
ever, such a correction factor may not be warranted if 
those sources detected at 160 /im are intr insically differ- 
ent from those not detected. Figure 18 shows the dis- 
tribution of 70 fim fluxes for these two subsamples. The 
sources detected at 160 fim are clearly weighted toward 



the bright end of the flux distribution. 

To further examine the properties of the sources not 
detected at 160 /im, we fit templates to their mean SED 
using a stacked 160 /im data point in several redshift and 
flux bins in the same manner as the detected sources. 
The resul ting comparison is shown as the colored points 
in Figure [19] While these stacked points show a similar 
amount ofscatter as the 160 /im detected sources, they 
do not show the same systematic offset. For the stacked 
sources, the luminosity obtained with the 160 /im point 
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Fig. 20. — Comparison of luminosities obtained with all of the available bands for the sources without 160 (im detections (left) and those 
with 160 (im detections (right) to the luminosities obtained from just the 24 (im data point. The dotted line represents the mean difference 
between the two and the cr value is the internal dispersion about that me an. The sources detected at 160 fim show a similar offset in their 
luminosities as they did in the 70/160 fim luminosity comparison (see Fig. |13| This is further indication that the sources detected at 160 fim 
are intrinsically different (i.e., have a colder dust component) than those that are not. 



is in agreement with the luminosity obtained without 
it. This suggests that the sources detected at 160 fim 
are indeed biased toward colder sources. Therefore, a 
correction to the Lir for the rest of the 70 /im sample 
is not necessary. Instead, we include the scatter in the 
relationship as an additional source of uncertainty and 
include this in our error bars in Figure |12| Thus the 
typical uncertainty for sources detected at 160 fim is 0.05 
dex while for those not detected it is ~ 0.2 dex. 

4.6. Estimates based on 24 fim 

Since MIPS is the most sensitive at 24 /im, samples 
selected at this wavelength are very useful and will con- 
tinue to be used in the future. Therefore, it is essential 
to understand how well their total infrared luminosities 
can be determined with just this data point. To test this, 
we obtained the best fit template using just the 24 fim 
data point. Since there is only one point, the templates 
are not really fit but instead the best match to the 24 /im 
flux at the source redshift is obtained from each of the 
template libraries without allowing for rcscali ng. For the 
final luminosit y estimate we chose to use the Siebenmor- 
gen & Kriigcl library since it was the best match for the 
vast majority of sources using 70 and 160 fim and has the 
largest variety of SED shapes. Figure [20] shows the com- 
parison between luminosity estimates using all available 
bands and using only the 24 /im point for both 160 fim 
detected and non-detected sources. The scatter in this 



relationship is much larger (cr ~ 0.4 dex) and the differ- 
ence can be up to 1 dex for some sources. A larger offset 
is also seen for the sources detected at 160 fim than those 
not, consistent with the results discussed in the previous 
section. The large scatter indicates that luminosity esti- 
mates for samples based solely on 24 fim detection alone 
can be off by a significant amount and many LIRG and 
ULIRG samples in the literature based on 24 fim may ac- 
tually be contaminated by lower luminosity sources since 
their number densities are higher though it is possible 
that this effect is mitigated by the slight offset (~ 0.1 
dex on average) obtained without the 70 fim data point. 

5. DISCUSSION 

In this paper we have presented a sample of 1503 70 fim 
sources in the COSMOS field and determined a reliable 
measure of their total infrared luminosity. We believe 
that this is the largest sample of 70 fim sources to date. 
When the source is detected at 160 fim we obtain an es- 
timate of Lir accurate to within 0.05 dex on average and 
0.2 dex for those without a 160 fim detection. This re- 
sult is much better than can be obtained from a 24 fim 
data point alone (~ 0.5 dex). In this section, we discuss 
in more detail the UV-to-FIR SEDs and the variations 
of the mean SED shape with Ltr. In particular we ex- 
plore three spectral regions (the UV, NIR, and MIR) 
since these have been found to be the most sensitive to 
the underlying power source (i.e., starbursts and AGN). 
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0.1 1.0 10.0 100.0 1000.0 

Rest-Frame Wavelength [/urn] 

Fig. 21. — Spectral energy distributions of all 1503 70 fim sources binned by luminosity and normalized at 1/im. The SEDs are color 
coded such that the ones with the bluest slope in the UV-optical are purple/blue and the reddest are red and the spectrum in between 
shows the range of SED colors. The median SED for each bin is overplotted in blue. Note the large variety in SED shapes, particularly 
in the UV, and the trends with Lir. At the highest luminosities there is an obvious trend toward power law SEDs all the way from the 
optical through the FIR. 
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Fig. 23. — Rest frame MIR spectral index (ctjjj) versus UV/optical spectral index (ag'l) divided into Ljr bins. Note that ci^q appears 
to span a larger range in values with increasing Ztr. 



We make use of the multiwavclcngth properties of the 
sample to identify a large population of AGN candidates. 

5.1. Analysis of SEDs 

Figures [6] and [21] show the full range of spectral en- 
ergy distributionsThat we observe for the 70 //m selected 
sample of galaxies binned by redshift and Ltr with the 
median SED in each bin over plotted as a blue line. 
The median SEDs are similar to sources observed lo- 
cally for sources at low luminosity (e.g., Dale etal. 2007: 
SINGS, iiR < lO n i ) and at high luminosity (e.g., 
U etal. 2009, in preparation: RBGS, L IR > 10 U L Q ). 
All of these studies have shown that the median UV-to- 
FIR SEDs have two prominent thermal peaks (the op- 
tical stellar peak and the FIR dust peak) and, in addi- 
tion, a range of spectral shapes, particularly in the UV 
(0.2-0.5 /mi), NIR (2-5 /im), and MIR (10-30 /im). With 
the large number of sources in our 70 /im sample we can 
finally begin to quantify the full extent of these varia- 
tions and explore how different parts of the SED corre- 
late with each other and with Ltr. As should already be 



known from studies of smal ler samples of local objects 
(e.g., RBGS: ISanders e t al.||2003| 1 Jy ULIRG Survey: 
|Kim fc Sanders||l998| , it is clear that a single archetyp- 
ical template for a given infrared luminosity (e.g., Arp 
220 or Mrk 231 for ULIRGs) would not be applicable for 
all objects with that luminosity. 

The median SEDs clearly show several general features. 
First, most of them show the 1.6 /im peak of the stellar 
bump. This feature becomes less pronounced at higher 
luminosities. The SEDs also display a prominent 5 /im 
dip at low luminosities which becomes less clear at high 
luminosities. The 7.7 /im PAH emission feature is also 
evident at Lir < 10 11 L Q near 8 /jm but becomes less ev- 
ident in the higher luminosity bins as rest frame 7.7 /im 
is no longer probed. However, in the highest luminos- 
ity bins, 24 /mi approaches rest-frame 7.7 /mi and this 
enhancement is still not observed, possibly due to the 
median SED being dominated by AGN at those luminosi- 
ties. Many of the SEDs show an increase in e mission to- 
ward the UV, indicative of the big blue bump dMalkan & 



(p v 

Sargent|1982| |Sanders et~aT|l989| [Elvis et al.||l994[ ) seen 
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Fig. 24. — Rest frame MIR spectral index (affj) versus NIR spectral index (c?|) divided into Ljr bins. Notice that at higher luminosities 
(I/rR > 10 11 Lq) a second branch appears with large values of a| (power-law SEDs) and intermediate values of ctfg- 



in many AGN. The GALEX data points, particularly at 
high redshift where they probe further into the FUV, 
point to a broad maximum for many of these sources. 

To provide a more quantitative measure of the spread 
in SED shapes, we compute three spectral indices (log- 
arithmic rest-frame vL v ratios) probing three different 
wavelength regimes: UV/optical from 0.2 to 0.5 /im 
(a#i), NIR from 2 to 5 /mi (a|), and MIR from 10 to 
30 fim (af[]). These spectral indices are plotted against 



each other in Figures |22|[24| binned by Ltr. All three of 



these plots show strong trends with Ltr. a\ increases 
with Ltr as sources become dominated by pow er-law 
SEDs (sec §5.2). This is also evident in Figure 24 where 
the sources with high a| form a separate distinct branch 
at intermediate values of For galaxies with low val- 
ues of all a w IS l° w a t l° w luminosities and steadily 
increases with L IR . a®^ spans a wide range at all L IR 
(a range of nearly 3 orders of magnitude) and has lower 
values on average at the most extreme luminosities. 



5.1.1. 160/70 /im Flux Ratios 

In addition to the SED shapes described above, our 
large 70 /im sample and the substantial number of de- 
tections at 160 /im makes it particularly useful to in- 
vestigate the 160/70 /im flux ratios as a crude mea- 
sure of the characteristic dust temperature for the FIR 
peak emission. Figure [25] shows the observed FIR colors 
(vF v (160) /vF v (70)) as a function of redshift and Ltr for 
all 70 /im sources detected at 160 /im. For those not de- 
tected at 160 /im, the colored points on the plot represent 
the colors obtained from their stacked 160 /im fluxes. To 
put these values in context we also over plot the flux ra- 
tios for local IRAS 60 /im selected sources from the RBGS 
over a wide range of L m (10 s - 10 12 5 i Q ) with MIPS 
observations obtained as a part of GOALS or SINGS 
(the Spitzer Infrared Nearby Gal axies Survey: Kenni- 
cutt et al.||2003| |Dale et al.||2007|) . A typical error bar 



tor the sources detected at 160 /zm is shown in the cor- 
ner. The overall scatter of the data points remains the 
same even with a cut at larger S/N, indicating that the 
observed trends are likely real. The local sources span 
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Fig. 25. — Observed far-infrared colors (vF 1/ (160)/vF 1 ,(70)) as a function of redshift (left) and Lin (right). A typical 1 a error bar is shown 
in the corner. The colored points represent the stacked data points (described in §4.4) for those sources not detected at 160 /im divided into 
several flux and redshift bins. Local infrared sources from the RBGS are over plotted with their FIR colors obtained from MIPS observations 
from SINGS (10 9 < L m < 10 11 Lq : open diamonds) and GOALS (10 11 < L m < 10 115 Lq: open triangle s; 10 115 < Lm < IP 1 25 Lp,: 
open squares). The gray lines represent the redshift evolution of empirical templates of local objects from Dale fc Helou| |2002| in the 
luminosity range 10 8 ' 3 < Lir < 10 13 ' 5 Lq. while the black lines highlight the luminosities at 10 9 , 10 10 , 10 11 , and 10 LZ Lq from top to 
bottom. 
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Fig. 26. — Same as Fig. |25| but converted into rest-frame colors using the best fit template to the infrared. Overall, the COSMOS sources 
span the same range in FIH colors as the local sources, however, at high luminosities there are a number of sources with colder SEDs than 
would be expected based on local sources. This may be evidence for a population of sources with colder SEDs at high redshift, however, at 
these redshifts the 160 ^m data point is not long enough to constrain the peak of the SED nor to obtain a precise dust temperature. The 
large r ang e in SED templates that fit these data points is indicated by the fact that these sources are the same ones with a large scatter 
in Fig. 



a large range in ^^(160) jvF v (70) from 0.2-5 and show 
a strong trend with Lir. Also ov er plotted in gray are 
a set of values obtained from the Dale & Helou (20021 
templates spanning an Ltr range of 10 S J — 10 iJ O from 
top to bottom as a function of redshift. These templates 
span most of the range occupied by the COSMOS 70 fim 
selected sources (the full set of libraries covers the entire 
range of colors observed at each redshift but we show 
just this one library for illustration). The sources not 



detected at 160 /im span the lower portion of the range 
occupied by the sources detected at 160 /im. This result 
is expected if the 160 /zm detection preferentially selects 
colder sources than the rest of the sample. 

We obtained rest-frame vF v (l&Q)/vF v (7Qi) values by 
using the 70 and 160 /im fluxes determined from the best- 
fit IR template used to estimate Ltr. These values are 
plotted in Figure [26] as a function of redshift and luminos- 
ity. As a note of caution, since these values are dependent 
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on the best-fit SED template this adds an extra source of 
uncertainty, particularly at the high redshift end (where 
shorter rest-frame wavelengths are observed). We find 
that the COSMOS sources span the same range of col- 
ors (and thus dust temperatures) as the local sources 
but that there is an apparent excess of colde r sources at 
higher lumin osities, similar to the result of |Symeonidis| 



et al. 



_ ( |2009| ). However, our interpretation is that this 
excess is likely a result of the longer wavelength selection 
(rest-frame 80 /jm at z ~ 1 versus 60 /im locally) and the 
fact that the peak of emission is not well constrained at 
these redshifts and thus a wide range of templates can fit 
the data points. Future surveys at FIR and submillimc- 
ter wavelengths with Herschel and SCUBA2 will be able 
to disentangle the temperature dependent selection ef- 
fects and investigate the dust temperature distributions 
of high redshift galaxies. 

5.2. AGN Diagnostics 

Both X-ray and radio detections are traditional meth- 
ods used to identify potential AGN in high-redshift 
galaxies. The vast majority of the 70 /im selected galaxies 
that are detected in the X-ray (86%) have an X-ray lumi- 
nosity greater than 10 42 ergs _1 and are therefore bona 
fide AGN. The 22 lower luminosity sources are all at low 
redshifts (z < 0.35). Most of the X-ray detected AGN 
have power law SEDs (though not all), particularly at 
high Lir . Many sources that are not detected as AGN in 
the X-ray also have a power-law SED. These are possibly 
heavily obscured (Compton t hick) AGN that are o paque 
even to hard X-ray emission ( |Polletta et a l. 2006). 

Since a\ is a measure of the v F v slope between 2 and 
5 fim it is therefore a good way to select power law galax- 
ies. Anything above zero on the y-axis of Figure [22] or 
the x-axis of Figure 24 has the shape of a power-law in 
the NIR, i.e., these objects do not have a strong 1.6 /mi 
bump nor a dip at 5 /im. The a\ > cutoff here corre- 
sponds to a F u spectral slope of 0.4 (similar to what is 
used as a power law slope cu toff in the literatu re, e.g., 
|Alonso-Herrero et al. 2006 and Donley et al. 2007 ) . X-ray 
detected AGIN span almost the entire area of the spectral 
index plots rather than being clearly separated the way 
the power law galaxies are. However, they do on average 
have lower afg as AGN in the local universe do. 

We select all sources with a? > to be power-law 
AGN candidates. A total of 166 (11.2%) objects satisfy 
this criterion. Forty-five of these are detected i n X-ray 
(~ 27%; a slightly low er fraction than fou nd by Alonso- 
Herrero et al.|2006| and |Donley et al.|2007 but consistent 



with the shallower COSMOS X-ray depths) . A few differ- 
ences between power-law galaxies detected in the X-ray 
and those not detected in the X-ray can be seen in their 
SEDs. Those sources detected in the X-ray tend to dis- 
play the full range of properties in the UV where only one 
or two sources that are not X-ray detected appear bright 
in the UV. The prominence of the UV emission in the X- 
ray detected power-law galaxies could be evidence that 
these sources tend to be less obscured than the non X-ray 
detected sources and hence are transparent to UV (and 
X-ray) emission. The SEDs of the non- X-ray detected 
power-law galaxies also tend to have a steeper slope. 

We use the infrared-radio correlation as another way 
to investigate the energy source. Using the limits defined 



in Yun et al. (2001 1 (5 times larger radio or IR flux than 
expected from the infrared-radio correlation) for radio 
excess (q < 1.64) and infrared excess (q > 3.0) objects we 
find many objects that fall into these two categories. In 
total there are 27 radio excess sources, 6 of which (22%) 
are also detected as X-ray AGN (as shown in Fig. 27 1. 
The fraction of objects that this population represents 
increases with L m from 4 (2%) at < 10 11 L , to 5 (2%) 
LIRGs and 16 (14%) ULIRGs. Radio excess sources are 
believed to be potential AGN hosts where the extra radio 
emission seen comes from either a compact radio core or 
from associated radio jets or lobes (Sanders & Mirabel 
1996). 

All 14 of the infrared excess objects are ULIRGs or 
HyLIRGs and represent 11% of the sample at these lu- 
minosities. 6 (43%) of them are X-ray detected AGN. 
A total of 25 of the obscured AGN candidates (64%) 
described in §5.2.2 are detected in the radio and have 
< q >= 2.54 ± 0.63. 5 of these (2 X-ray AGN) are in- 
frared excess objects. It is interesting to note that the IR 
excess objects contain a higher fraction of X-ray detected 
AGN than the radio excess objects. It is also interest- 
ing to note that many of the radio and infr are d excess 
sources also have a power-law SED (see Fig. 27 1. 



5.2.1. IRAC Color Selected AGN 
Color selection in the mid-infrared has also been shown 



to be an effective way to sel ect potential AGN (Stern 



et al. 2005 Lacy et al. |2004[ ) as these colors can sepa- 



rate star-forming galaxies with a stellar bump at 1.6 /im 
from those with a power-law like SED. Also, this method 
can potentially find dust obscured AGN since the MIR 
wavelengths are not as affected by ex tinction as the UV 
and optical. Following the criteria of |Stern et al. (2005) 
we used the color-color IRAC plot (see Fig. [28| to se- 
lect potential AGN and compare these sources with those 
found using four other methods. A total of 232 sources 
meet the Stern et al. criteria, representing ~ 15% of the 
entire 70 /jm sample. Of these 232, 82 (35%) are de- 
tected as AGN in the X-ray. This is 62% of the total 
X-ray AGN population, which means the other 50 X- 
ray AGN are not selected by this method. According 
to Stern et al. , this method is good at selecting broad 
line AGN (selects ~ 91%) and not as good at selecting 
narrow line AGN ( ~ 40%) with an overall contamina- 



tion rate of ~ 20%. Donley et al. (20071 show that this 



region of the plot can be contaminated by starburst dom- 
inated ULIRGs (like Arp 220) and normal star-forming 
galaxies. If we assume that 20% of the 232 objects are 
potentially contaminants, this leaves us with 186 actual 
AGN, or about 1.5 times the number detected in the X- 
ray. Of the potentially obscured objects described in the 
next section (F(24//m)/F(.R) > 10 3 ), 36 out of the 60 
match these color criteria as do 33 out of the 56 that also 
have R — K> 4.5. Most of t he power law galax ies (62%) 



fall in a narrow region in the Stern et alT] ( |2005[ ) selection 
box but many (~ 35%) fall outside of tnis region. Very 
few of these objects fall in the lower left corner of the se- 
lection region, indicating that this is possibly the region 
that contributes most of the contamination in th e selec - 

[20081). 



tion (consistent with the results of Donley et al. 
13 of the radio excess sources are also selected 
Stern et al. criteria. 



jy the 
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Fig. 27. — Radio/infrared flux ratio (q) vs. Lir (same as Figure 1 10 [ ) with X-ray detected AGN (red), objects with power-law SEDs 
(blue), IRAC color selected AGN candidates (magenta), and obscured AGN candidates (cyan) over plotted from top to bottom. The solid 
line indicates the median value of q (2.34) for the entire sample and the dotted lines represent the dividing lines for radio (q < 1.64) and 
infrared (q > 3.0) excess sources following Yun etal. (2001). The panels on the right show the distribution of the entire sample with radio 
detections and the different AGN selections over plotted in color. 
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([5.8] - [8.0]) ve ,„ 

Fig. 28. — |Stern et al.| | |2005| l IRAC color-color plot for selecting AGN candidateswith X-ray detected AGN (red), obj ects with power-law 
SEDs (blue), radio excess objects (green), and obscured AGN candidates (cyan) over plotted from top to bottom. The |Stern et al.] ( |2005[ l 
selection box is represented by the dashed lines. 
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Fig. 29. — MIR (24 fim) to optical (R-band) flux ratio as a function of R — K color with X-ray detected AGN (red), objects with 
power-law SEDs (blue), radio excess objects (green), and IRAC color selected AGN candidates (magenta) over plotted from top to bottom. 
The dotted lines represent the F(24)/F(R) > 10 3 and R — K > 4.5 selection cuts from Fiore etal. 2008 for obscured AGN candidates. 
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5.2.2. Obscured AGN Candidates 



Fiore etal. (2008 



|2009[ ) have shown that obscured 
AGN tend to have high MIR to optical flux ratios and 
red R — K colors. The 24 /im to R-band flux ratio for 
each of the 70 /im sources in the sample is shown in 
Figure [29] as a function of R — K color with AGN se- 
lected through four other methods over plotted. The 
dotted lines show the cuts of ( R — K)y eKei > 4.5 and 



Fiore et al. (2008k for 



at high redshift . 56 of 



F(24/xm)/F(R) > 10 3 defined in 
selecting candidate obscured AG1S " 
the galaxies in our sample make both of these cuts and 14 
of them are detected in the X-rays. Three are detected 
only in the hard band and one only in the soft band. 
The remaining sources have a median hardness ratio of 
0.2. The redshifts of these galaxies range from 0.7 to 2.5 
with the median b eing z = 1.5, confirming the results of 



Fiore et ah] ( 2008 1. Additionally, four galaxies match the 
MIR to optical flux ratio criteria but are bluer than the 
(R — -ftT)vcga > 4.5 magnitude cutoff. Of these, one is 
detected as an AGN in the X-ray with a detection only 
in the hard band. These four sources have redshifts that 
range from z = 1.2 to 2.4 with a median redshift of 2.1 - 
a similar distribution to the sources with red R — K col- 
ors. The 24 /im to R-band flux ratio correlates with the 
total infrared luminosity (as expected for obscured AGN 
where the luminosity from the nuclear regions is greatly 
reduced by the obscuring material) and those above the 
cutoff of 10 3 are predominantly ULIRGs and HyLIRGs. 
We consider all 61 of these galaxies to be obscured AGN 
candidates. 

It is also worth noting that this selection is similar to 
other IR/optical color sel ections in the lit erature, such 



bright, optically faint galaxies 1 


Yan et al. 


2004 


Houck 


et al. 2005 Weedman et al. 


20( 


16 


). The d 


ust obscured 



galaxy color cut {R — [24J > 
the one used here and is expected to select high redshift 
[z ~ 2) extremely luminous galaxies with large column 
densities of dust. 

5.2.3. AGN Fraction Among Luminous and Ultralummous 
Infrared Galaxies 

In the local universe, the AGN fraction of galaxies is 
known to increase very strongly with Ltr. Using spectro- 
scopic diagnostics o f a sample of 200 lu minous infrared 
galaxies from IRAS, |Veilleux et al.l(|1995| ) found that 62% 
of objects at the highest luminosities wer e AGN. A sim- 
ilar s tudy of the 1 Jy sample of ULIRGs ( Veilleux et al. 
1999] fou nd that above Lm > 10 12 3 ~ 50% of objects 
are AGN. Tran et al. (2001 1 found a similar transition lu- 
minosity of -Lir > lO 1 ^ 4 Lq between objects whose lumi- 
nosities are dominated by star-burst emission and those 
dominated by AGN. Above this transition they find that 
"most ULIRGs appear AGN- like" . With the large sam- 
ple of objects over the wide range of luminosities probed 
by our 70 /im selected sample, we can address the ques- 
tion of whether this trend holds out to higher redshifts. 

Table [5] shows the total number of galaxies in each Ltr 
bin along with the number of X-ray detected sources over 
the full area, in the area covered by Chandra, and in the 
area with the deepest Chandra coverage. The total frac- 
tion of X-ray detected AGN increases with increasing 
L m , from 5% at L m < 10 10 L Q to - 15% for ULIRGs 




log(L«) 

Fig. 30. — Fraction of 70 /im sources selected as AGN candidates 
by a ny one of the five different selection methods illustrated in 
Figs. [27j|29| as given in Column 9 of Table [6] 1 a Poissonian error 
bars are shown. Many of these sources are selected by more than 
one method though each source is counted only once. The fraction 
goes from < 5% at L m < 10 11 L to 100% at Ltr > 10 13 Lq. 



and ~ 42% for HyLIRGs, however, these numbers are 
lower limits due to the uneven X-ray coverage across the 
field. The fraction increases when limited to the smaller 
areas with deeper X-ray coverage, though the small num- 
bers make a statistical analysis difficult. This indicates 
that as many as 12% of LIRGs, 30% of ULIRGs, and 43% 
of HyLIRGs are AGN with detectable X-ray emission. 
For the rest of the discussion here, we use the fractions 
obtained over the entire area. 

Table [6] shows the fraction of candidate AGN found by 
each method (including X-ray detected, radio selected, 
IRAC color selected, power-law galaxies, and potentially 
obscured objects) . If we add up all of the galaxies that 
fall into any of the possible AGN candidate categories 
(total column in Table |6| then we find that the AGN 
fraction ranges from 3% at £ir < 10 10 Lq to ~ 70 — 



100% for ULIRGs and HyLIRGs (illustrated in Fig. (30 1. 
There is likely to be some contamination in this estimate, 
particularly from the IRAC color selection method. If we 
leave out the IRAC selected objects, the AGN fraction 
goes from 2% for objects with Ltr < 10 10 Lq to 10% for 
LIRGs, 51% for ULIRGs, and 97% for HyLIRGs. This 
result is in excellent agreement with the AGN fraction in 
the local universe, indicating that high redshift ULIRGs 
are indeed similar to local sources in this respect. 

The redshift and Ltr distribution for the AGN candi- 
dates found using each method is shown in Figure [31] 
Those objects that are also selected in the X-ray are 
shown by the filled histogram. It is interesting to note 
that while the X-ray selection finds objects at a wide 
range in redshifts and Ltr, each of the other methods 
has a narrower range and typically identifies objects at 
higher luminosities. The IRAC color selection finds two 
peaks in the redshift distribution, one at z ~ 0.35 and the 
other at z ~ 1.15. For the radio, IRAC, and power-law 
selections the X-ray seems to detect objects over the full 
range of redshifts and Ltr. The majority of the power- 
law objects are ULIRGs or HyLIRGs with only a few 
(~ 10) with LIRG luminosities. 
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Fig. 31. — Redshift (left) and infrared luminosity (right) histograms for AGN selected by each of the five different selection methods: 
(from top down) X-ray selected, radio excess sources, IRAC color selected sources, sources with power law SEDs, and obscured AGN 
candidates. The filled histograms represent the objects in each method that are also detected in the X-ray. The median Lir is plotted for 
each method. 
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6. SUMMARY 

We have presented a large robust sample of 1503 un- 
confused 70 /im selected galaxies with S/N > 3 (ap- 
proximate flux limit of 6.5 mJy) from the 1.8 deg 2 ACS- 
COSMOS field. These sources span a redshift range of 
0.01 < z < 3.5 with a median redshift of 0.5 and an in- 
frared luminosity range of 10 8 < Lm < 10 14 L© with a 
median luminosity of 10 11 ' 4 L© This is the largest sam- 
ple of 70 /xm sources to date. The large multiwavelength 
data set available for the full COSMOS field has allowed 
us to compile the complete spectral energy distributions 
for each object from the ultraviolet to the far-infrared 
and obtain an estimate of the total infrared luminosity. 
Spectroscopic redshifts are available for ~ 40% of the 
sample and for the remainder we use photometric red- 
shifts with a precision of 0.02 x (1 + z). The coverage 
in the X-ray, radio, and IRAC bands have allowed us to 
identify a large sample of potential AGN. From the anal- 
ysis of the multiwavelength properties of these objects, 
the following conclusions can be drawn: 

1. The long wavelength Spitzer-MIPS 70 and 160 /im 
data are essential for obtaining a more accurate 
and reliable estimate of the total infrared lumi- 
nosity than was previously possible using 24 /im 
data alone. The typical uncertainty in Lm when 
160 /im is detected is ~ 0.05 dex while when it is 
not detected it is ~ 0.2 dex. Using only the 24 /im 
data point we find the dispersion is significantly 
increased to ~ 0.5 dex. 

2. The shape of the mean SED for the 70 /im selected 
sample including 18 different bands from the FUV 
through the FIR is similar to what has been ob- 
served for local LIRGs and ULIRGs. In particular, 
the maximum IR to optical ratio (~ 100) and the 
maximum Lm (~ 10 14 L@) is the same as seen lo- 
cally. Our large sample presents a more complete 
picture of the range of SED shapes, particularly in 
the UV where the ratio of the optical (0.5 /im) to 
ultraviolet (0.2 /im) luminosities varies by nearly a 
factor of 1000. Similarly, the ratio of luminosities 
in the near (2 to 5 /im) and mid-infrared (10 to 
30 /im) varies by more than a factor of 100. Al- 
though our selection at 70 /im is biased toward a 
warmer sample at higher redshift it still appears 
that the dispersion of shapes is similar at all red- 
shifts. Clearly, one needs to consider a wider range 
of spectral templates as opposed to the typical cold 
(e.g., Arp 220) and warm (e.g., Mrk 231) templates 
used to represent infrared galaxies. 

3. A comparison of our luminosity estimates using 
model fits with and without the 160 /im data points 
indicates that the sources detected at 160 /im have 
a colder SED than those not detected. This is con- 
firmed by the difference in the FIR colors between 
the two subsamples. A comparison with local in- 
frared selected sources finds the same range in FIR 
colors overall but also presents evidence that high 
redshift and high luminosity sources detected at 
160 /im in our sample are colder than local sources 
of comparable luminosity. This is likely a combi- 
nation of a selection effect due to the longer wave- 



length 160 /im detection and the poor constraints 
that 160 /im places on the peak of emission and 
dust temperature at these redshifts (z > 1.5). Fu- 
ture FIR and submillimctcr studies with Herschel 
and SCUBA2 will be necessary to disentangle these 
possible causes. 

4. We use five of the most commonly used multi- 
wavelength selection techniques to identify poten- 
tial AGN in our sample, including 132 X-ray de- 
tected AGN, 27 sources with radio fluxes five times 
greater than that expected from the infrared-radio 
correlation, 166 sources with power-law SEDs, 232 
IRAC color selected galaxies, and 61 potentially 
obscured (F(24//m)/F(R) > 10 3 ) AGN. We find 
that there is significant overlap among these differ- 
ent methods resulting in a total of 354 sources with 
signatures of an AGN. 

5. We find that the total fraction of AGN increases 
strongly with Ljr similarly to the behavior of 
galaxies in the local universe. Our estimate indi- 
cates that nearly 70% of 70 /im selected ULIRGs 
and all HyLIRGs contain a powerful AGN. Less 
than half of the candidate AGN selected based on 
several different methods are detected in the X-ray 
providing evidence for a significant population of 
obscured AGN among this infrared selected sam- 
ple. Since a 70 /im selection selects warmer SEDs 
at higher redshifts it remains to be seen whether 
this AGN fraction is representative of all ULIRGs 
at these redshifts. 

The complete 70 /im selected catalog of all 1503 sources 
including their coordinates, redshifts, and total infrared 
luminosities (Table J3]) will be made publicly available 
with the online version of this paper. 
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TABLE 1 

Summary of COSMOS Multiwavelength Data Sets 



Instrument 


Band 


Date of Observations 


Depth (5<r) 


# in 70 (im sample 


opirzer/ ix\./\v> 


3.6 /im 


zuuo 


n Q 1 1 Tir 

u.y f-ioy 


10UO 


Spitzer/IRAC 


4.5 //m 


2006 


1.7 ujy 


1503 


Spitzer /IRAC 


5.6 |ini 


2006 


11 3 Tv 

1 J- »V 1AU V 


1501 


Spitzer/IRAC 


8.0 /im 


2006 


14.6 /ijy 


1503 


Spitzer/MIPS 


24 /im 


2006-2008 


60 /ijy 


1503 


Spitzer/MIPS 


70 fim 


2006-2008 


4.5 mjy 


1503 


Spitzer/MIPS 


160 /im 


2006-2008 


65 mjy 


463 


HST/ACS 


F814W 


2003-2005 


27.2 mag 


1503 


HST/NICMOS 


F160W 


2003-2005 


25.9 mag 


74 


Galex 


FUV.NUV 


2004 


25.7, 26.0 mag 


892, 1284 


XMM 


0.5-10 keV 


2004-2005 


5 X 10 -16 erg cm -2 s _1 


119 


Chandra 


0.5-10 keV 


2006-2007 


5.7 X 10 -16 erg cm -2 s _1 


103 


CFHT/MEGACAM 


u* 


2003-2007 


26.5 mag 


1502 


Subaru/SuprimeCam 


BjVjg+r+i+z+ 


2004-2005 


j + ~ 26.2 mag 


1503 


UKIRT / WFC AM 


J 


2004-2007 


23.7 mag 


1498 


CFHT / WIRC AM 


K S 


2005-2007 


23.5 mag 


1503 


VLA 


1.4 GHz 


2003-2005 


55 /ijy 


562 



TABLE 2 

Spectroscopic Redshift Surveys 



Survey Reference Instrument Wavelength Range [A] # in 70 (im sample 



24 (im sources 


Kartaltepe et al. , in prep. 


Keck II/DEIMOS 


4000-10000 


62 




Capak, Salvato et al. , in prep 


Keck II/DEIMOS 


4000-10000 


18 


zCOSMOS 


Lilly et al. 2009 


VLT/VIMOS 


5550-9650 


355 


XMM sources 


Trump et al. 2007, 2008 


Magellan/IMACS 


5600-9200 


223 


Quasar candidates 


|Prescott et al.|2006| 
TlnEuita et al. 2UL)y 


MMT/Hectospec 


3100-9000 


43 




VLT/FORS1 


3300-11000 


4 


SDSS 


|Abazajian et al.|2009| 


SDSS Spectrograph 


3900-9100 


65 


2dFGRS 


|Oolless et al.|2U01| 


2dF Spectrograph 


3500-10000 


34 



TABLE 3 

COSMOS 70 jitM Selected Source Catalog 



ID 


Name 


RA 


Dec 


24 /im flux a 


70 (an flux 160 /im flux 


Redshift 


Flag b 


log(L IR /L ) c 


1017 


COSMOS J100046.07+013439.89 


150.19197 


1.57775 


0.69 + 0.02 


7.25 + 2.13 


0.518 


I 


11.16 + 0.25 


1077 


COSMOS J100034.22+013543.18 


150.14259 


1.59533 


1.08 + 0.29 


16.03 + 2.03 


0.23 


P 


10.87 + 0.22 


1081 


COSMOS J100133.37+013547.59 


150.38905 


1.59655 


0.75 + 0.04 


10.82 + 2.00 


0.23 


P 


10.44 + 0.22 


1082 


COSMOS J100041. 55+013552. 68 


150.17314 


1.59797 


1.62 + 0.35 


22.26 + 2.45 


0.37 


P 


11.26 + 0.23 


1083 


COSMOS J100241. 78+013550. 07 


150.15190 


1.59699 


0.14 + 0.02 


18.61 + 2.12 


1.20 


P 


12.62 + 0.32 



a All fluxes and errors are in mjy. Fluxes and errors for all of the bands discussed in this paper arc presented in the online version of the paper. Here, only 
the MIPS data points are shown as an example. 

b Redshift Flag — D: DEIMOS spectroscopy on Keck II, Z: zCOSMOS spectroscopy from VLT/VIMOS, I: IMACS spectroscopy from Magellan, S: SDSS 
spectroscopy, 2dF: 2dFGRS spectroscopy, F: FORS1 spectroscopy, M: Hcctospcc spectroscopy from MMT, P: Photometric redshift from Ilbcrt et al. 2009 or 
Salvato et al. 2009 

c Total infrared luminosity and uncertainties estimated from the template fits described in §4. 



70 Sources: SEDs and Lir 



TABLE 4 

Offset in L ir Estimate Without 160 /xm Data Point 

Flux bin [mjy] z bin med(fiux) [mjy] med(z) # of sources Mean Offset Dispersion 

0.4-0.8 ■■■ 6.6 0.84 61 0.41 0.20 

0.8-12.0 •■■ 9.4 0.58 117 0.33 0.16 

12.0-20.0 ■■■ 14.7 0.41 140 0.20 0.18 

> 20.0 ••• 32.8 0.22 141 0.09 0.14 



0.0- 0.2 24.8 0.12 106 0.17 0.16 

0.2 - 0.4 15.9 0.27 124 0.14 0.15 

0.4- 0.6 12.8 0.49 69 0.19 0.22 

0.6- 0.8 12.2 0.69 49 0.28 0.20 

0.8- 1.0 11.4 0.90 50 0.30 0.19 

1.0- 1.5 10.6 1.15 39 0.31 0.25 

1.5 - 3.0 8.3 1.78 22 0.32 0.24 



TABLE 5 

Fraction of 70 fiM Sources Detected in the X-ray 



log(£m/£o) 


Median 


Number 


X-ray: Full area a 


X-ray: Chandra area b 


X-ray: Chandra deep area c 




Redshift 




# (%) 


# (%) 


# (%) 


< 9.0 


0.02 


3 











9.0 - 10.0 


0.10 


58 


3 (5.2) 


1 (3.6) 





10.0- 10.5 


0.17 


140 


8 (5.7) 


4 (5.2) 


2 (4.9) 


10.5 - 11.0 


0.27 


281 


18 (6.4) 


14 (8.0) 


9 (9.6) 


11.0 - 11.5 


0.47 


376 


35 (9.3) 


27 (12.0) 


15 (13.6) 


11.5 - 12.0 


0.77 


311 


28 (9.0) 


27 (14.1) 


14 (14.7) 


12.0 - 12.5 


1.05 


219 


29 (13.2) 


19 (13.0) 


14 (16.3) 


12.5 - 13.0 


1.47 


84 


20 (23.8) 


13 (29.5) 


9 (34.6) 


> 13.0 


2.00 


31 


13 (41.9 ) 


6 (46.2) 


3 (42.9) 


All Lir 


0.50 


1503 


154 (10.2) 


111 (12.3) 


66 (14.0) 



a Percentage of 70 (im sources detected in the X-ray by XMM and Chandra over the full ACS area of 
the COSMOS field. The XMM observations cover 100% of this area. 

b Percentage of 70 /im sources detected in the X-ray by XMM and Chandra in the area of the field 
observed by Chandra (55% of the ACS area) as illustrated in Figure |7| 

c Percentage of 70 /im sources detected in the X-ray by XMM and Chandra in the area of the field with 
deep Chandra coverage (30% of the ACS area). 



TABLE 6 
AGN Fraction by Luminosity Bin 



log(L m /L ) 


Median 


Number 


X-ray AGN a 


PL b 


IRAC C 


Obscured d 


Radio Excess 


Total* 


Redshift 




# (%) 


# (%) 


# (%) 


Candidates 


# (%) 


# (%) 


< 9 


0.02 


3 








1 (33) 








1 (33) 


9.0 - 10.0 


0.10 


58 








1 (1.7) 





2 (3.4) 


2 (3.4) 


10.0 - 10.5 


0.17 


140 


2 (1.4) 





2 (1.4) 





1 (0.7) 


5 (3.6) 


10.5 - 11.0 


0.27 


281 


12 (4.3) 





10 (3.6) 





1 (0.4) 


18 (6.4) 


11.0 - 11.5 


0.47 


376 


28 (7.4) 


3 (0.8) 


19 (5.1) 


1 (0.3) 


3 (0.8) 


40 (10.6) 


11.5 - 12.0 


0.77 


311 


28 (9.0) 


8 (2.6) 


40 (12.9) 


9 (2.9) 


2 (0.6) 


61 (19.6) 


12.0 - 12.5 


1.05 


219 


29 (13.2) 


73 (33.3) 


82 (37.4) 


22 (10.0) 


11 (5.0) 


123 (56.2) 


12.5 - 13.0 


1.47 


84 


20 (23.8) 


52 (61.9) 


56 (66.7) 


18 (21.4) 


5 (6.0) 


73 (86.9) 


> 13.0 


2.00 


31 


13 (41.9) 


30 (96.8) 


21 (67.7) 


11 (35.5) 


2 (6.5) 


31 (100.0) 


All Lir 


0.50 


1503 


132 (8.8) 


166 (11.0) 


232 (15.4) 


61 (4.1)) 


27 (1.8) 


354 (23.6) 



a Number and percentage of objects that are detected as AGN in the X-ray (see §6.4 for details). 
b Number and percentage of objects that have power-law SEDs as defined in §6.4 



(2005) 



match the 



c Number and percentage of objects whose colors match the AGN candidate criteria of Stern et al. 
d Number and percentage of objects whose infrared to optical flux ratios (E(24 fim) /F(R) > 
Fiore et al. 2008 criteria of potentially obscured AGN. 

c Number and percentage of objects that with a radio excess (5 times more radio flux than expected from the 
radio-IR correlation). 

f Total number and percentage of objects that fall into any of the AGN candidate categories. 



